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I. INTRODUCTION 



It is a well known historical fact that the 1925-1926 discoveries of two equivalent 
mathematical formulations of quantum mechanics — Heisenberg's matrix form followed by 
Schrodinger's wave mechanical form — preceded the development of a physical interpretation 
of these formalisms jl|. The first important ingredient of the conventional interpretation was 
Born's 1926 identification of the squared modulus of a complex Schrodinger wavefunction 
as a probability Q. The second ingred.ent developed in !927 was Heisenberg's Uncertainty 
Principle (UP) [3]. To these may be added Bohr's Complementarity Principle which has a 
more philosophical flavour (4 ] . 

Heisenberg's original derivation of his position-momentum UP combined the formula for 
the resolving power of an optical microscope extrapolated to a hypothetical gamma ray 
microscope, with the energy and momentum relations for a single photon, in analysing the 
inherent limitations in simultaneous determinations of the position and momentum of an 
electron. His result indicated the limits of applicability of classical notions, in particular the 
spatial orbit of a point particle, in quantum mechanics. 

More formal mathematical derivations of the UP, using the Born probability interpreta- 
tion, soon followed. Prominent among them are the treatments of Kennard, Schrodinger, 
and Robertson fl. Such a derivation was a.so presented by Heisenberg in his 19 30 Chicago 
lectures [6]. 

The Heisenberg position- momentum UP is basically kinematical in nature. In contrast, 



the Bohr UP for time and energy involves quantum dynamics in an essential manner 
Later work on the UP has introduced a wide variety of ideas jq] and interpretations of the 



fluctuations or the uncertainties involved [9], such as in entropic 10J and other formula- 
tions 

Even for a one-dimensional quantum system, the Schrodinger- Robertson form of the UP 
displays more invariance than the Heisenberg form. Thus while the latter is invariant only 
under reciprocal scalings of position and momentum, and their interchange amounting to 
Fourier transformation, the former is invariant under the three-parameter Lie group Sp(2, R) 
of linear canonical transformations. Fourier transformation, as well as reciprocal scalings, 
belong to Sp(2, R) Q. The generalisation of the Schrodinger-Robertson UP to any finite 



number, n, of degrees of freedom displays invariance under the group Sp(2n, R) [13J 



The purpose of this paper is to outline an invariant theoretic approach to general un- 
certainty relations for quantum systems. It combines a recapitulation and reexpression of 
some past results with some new ones geared to practical applications. The analysis 
throughout is in the spirit of the Schrodinger-Robertson treatment, and, in particular, our 
considerations do not cover the entropic type uncertainty relations. All our considerations 
will be kinematical in nature. 

The material of this paper is presented as follows. Section II sets up a general framework 
and procedure for deriving consequences of the positive semidefmiteness of the density matrix 
of a general quantum state, for the expectation values and fluctuations of a chosen (linearly 
independent) set of observables for the system. This has the form of a general uncertainty 
relation. A natural way to separate the expressions entering it into a symmetric fluctua- 
tion part, and an antisymmetric part contributed by commutators among the observables, 
hence specifically quantum in origin, is described. With respect to any unitary symmetry 
operation associated with the system, under which the chosen observables transform in a 
suitable manner, the uncertainty relation is shown to transform covariantly and to be pre- 
served in content. In Section III this general framework is applied to the case of a quantum 
system involving n Cartesian canonical Heisenberg pairs, i.e., an n-mode system ; and to the 
fluctuations in canonical 'coordinates' and 'momenta' in any state. The resulting n-mode 
generalization of the original Schrodinger-Robertson UP is seen to be explicitly covariant 
under the group Sp(2n, R) of linear homogeneous canonical transformations. Section IV 
returns to the single mode system, but considers as the system of observables the infinite set 
of operator polynomials of all orders in the two canonical variables. The treatment is formal 
to the extent that unbounded operators are involved. An important role is played by the set 
of all finite-dimensional real nonunitary irreducible representations of the covariance group 
Sp(2, R). We follow in spirit the structure of the basic theorems in the classical theory 
of moments. Thus the formal infinite-dimensional matrix uncertainty relation is reduced 
to a nested sequence of finite-dimensional requirements, of steadily increasing dimensions. 
While this case has been treated elsewhere 141 ]. some of the subtler aspects are now carefully 



brought out. In this and the subsequent Sections the method of Wigner distributions is used 
as an extremely convenient technical tool. Section V treats in more detail the uncertainty 
relations of Section IV that go one step beyond the original Schrodinger-Robertson UP. Here 
all the fourth order moments of the canonical variables in a general state are involved. Their 



fully covariant treatment brings in the denning and some other low dimensional represen- 
tations of the three-dimensional Lorentz group 5*0(2, 1). It is shown that the uncertainty 
relations (to the concerned order) are all expressible in terms of SO (2, 1) invariants. In 
Section VI we describe an interesting aspect of the Schrodinger-Robertson UP in the light of 
three-dimensional Lorentz geometry, which becomes particularly apparent through the use 
of Wigner distribution methods. We argue that this should generalise to the conditions on 
fourth (and higher) order moments as well. The paper ends with some concluding remarks 
in Section VII. 



II. GENERAL FRAMEWORK 



We consider a quantum system with associated Hilbert space "H, state vectors 
■ • • and inner product ((f)\ip) as usual. A general (mixed) state is determined by a density 
operator or density matrix p acting on "H and obeying 



p>0, Trp=l. 



(2.1) 



Then Tr p 2 = 1 or < 1 distinguishes between pure and mixed states. Any hermitian observ- 
able A of the system possesses the expectation value 



(A)=Tr(pA) 



(2.2) 



in the state p, the dependence of the left hand side on p being generally left implicit. 

We now set up a general method which allows the drawing out of the consequences of the 
nonnegativity of p in a systematic manner. This along with two elementary lemmas will be 
the basis of our considerations. 

Let A a , a — 1, 2, • • • , N be a set of TV linearly independent hermitian operators, each 
representing some observable of the system. We set up two formal A-component and (A+l)- 
component column vectors with hermitian operator entries as follows : 

I i \ 



A 



A 




\ An J 



A l 



(2.3) 



\ An J 



From A we construct a square (N + l)-dimensional 'matrix' with operator entries as 



/ 



n = aa 7 



1 



A, 



A a A b 



\ 



(2.4) 



/ 



Since (A a A b )^ = A b A a , Cl is 'hermitian' in the following sense : taking the operator hermitian 
conjugate of each element and then transposing the rows and columns leaves Cl unchanged. 
In a state p we then have an (N + l)-dimensional numerical hermitian matrix Q of the 
expectation values of the elements of Cl : 



Q= (Q) = Tr(pfi) = 
i.e., tt ab = Ti(pCl ab ) ; 

fit = a 



( 1 

V : 



( A b) 



(A a A b ) 



\ 



(2.5) 



Now for any complex (N + 1) component column vector C = (cq, c±, ■ ■ ■ , Cjv) t we have 

c j nc = cU(cU) f > 0, 

(C^C) = C f fiC > 0, (2.6) 

leading immediately to : 

Theorem 1 Positivity of p imputes positivity to the matrix Q, for every choice of A: 



p > Q = (Cl) = Tr(ptt) > 0, V^. 



(2.7) 



This is thus an uncertainty relation valid in every physical state p. 

Remark : It is for the sake of definiteness and keeping in view the ensuing applications that 
we have assumed the entries A a of A and A to be all hermitian. This can be relaxed and 
each A a can be a general linear operator pertinent to the system. The only change would 
be the replacement of A T in Eq. (\2A\i by ./v, leading to a result similar to Theorem 1. 

Depending on the basic kinematics of the system we can imagine various choices of the A a 
geared to exhibiting corresponding symmetries or covariance properties of the uncertainty 



relation (12.7ft . Specifically suppose there is a unitary operator U on H such that under 
conjugation the A a go into (necessarily real) linear combinations of themselves : 

uu ] = U ] U = 1, 
U' 1 A a U = R ab A b , 
U~ 1 AU = KA, 

(i o\ 

re = , re = re, 6 . (2.8) 
\ 0R J 

The matrix i? here is real iV-dimensional nonsingular. Then combined with Eq. 02.5 j) we 
have : 

p' = UpU~ 1 =>Q' = Tr(p'Q) = Tr{pU~ l AA T U) 

= kqk t , 

Q>0^Q'>0. (2.9) 

This is because the passage Q — > fi' is a congruence transformation. Thus the uncertainty 
relation (12 ,7p is covariant or explicitly preserved under the conjugation of the state p by the 
unitary transformation U. 

We now introduce two lemmas concerning (finite-dimensional) nonnegative matrices, 
whose proofs are elementary : 

Lemma 1 For a hermitian positive definite matrix in block form, 

Q = Q ] = | A ° | , ^ = A , £t = B, (2.10) 

\c B) 

we have 

Q >0 A> and B -C A~ l C ] > 0. (2.11) 

The proof consists in noting that by a congruence we can pass from Q to a block diagonal 
form [ijj] : 

1 \ I A W 1 ', 

Q=\ \ . (2.12) 

-CA- 1 1 / \ B-CA- 1 ^ I \ -CA- 1 1 
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Lemma 2 If we separate a hermitian matrix Q into real symmetric and pure imaginary 
antisymmetric parts R, S then 

Q = gt = R + iS > o, det,S ^ R > 0. (2.13) 

The nonsingularity of S means that Q must be even dimensional. (The proof, which is 
elementary, is omitted). 

Now we apply Lemma 1 to the (N + l)-dimensional matrix Q in Eq. (12.51) . choosing a 
partitioning where B is iV x N, C is N x 1 and is 1 x N : 

/ A Ct \ 

n = : A = 1, B = ((A a A b )), C = {{A a )). (2.14) 
\C B) 

Then from Eq. (12.111) we conclude : 
Theorem 2 

Q = (((Aa - (A a ))(A b - (A b )))) > 0. (2.15) 

All expectation values involved in the elements of the N x N matrix Cl are with respect to 
the state p. 

The motivation for the definitions of A, Cl as in Eq. (12.31) is now clear : after an application 
of Lemma 1 we immediately descend from the matrix Q to the matrix Cl involving only 
expectation values of products of deviations from means. It is then natural to write the 
elements of Cl as follows : 

AA a = A a - (Aa), 

Cl ab = (AA a AA b ). (2.16) 

We revert to this form shortly. 

The covariance of the statement (12.151) . Theorem 2, under a unitary symmetry U acting 
as in Eq. ( 12.8ft follows from a brief calculation : 

p -> p' = UpU 1 =^ 
U'\A a - Tr(p'A a ))U = U~ l A a U - ^(pU^AJJ) 

= R ab (A b -Tr(pA b )); 
Cl^Cl' = RClR T ; 
Cl>0^Cl'>0. (2.17) 



We now return to Eq. (I2.16p . The state p being kept fixed, we can split the hermitian 
N x N matrix Q into real symmetric and pure imaginary antisymmetric parts as follows : 

^ab = V ab (p ; A) + - u ab (p ; A), 

V ab {p- A) = V ba {p- A) = ^({AA a , AA b }) 

= \{{A a , A b })-(A a )(A b ); 
u ab (p ; A) = -io ba (p ; A) = -i ([A a , A b ]}. (2.18) 

The brackets [ • , • ] and { • , • } denote, as usual, the commutator and anticommutator 
respectively. The natural physical identification of the N x N real symmetric matrix 
V(p; A) = (V ab (p; A)) is that it is the variance matrix (or matrix of covariances) asso- 
ciated with the set {A a } in the state p. The uncertainty relation (I2.15P now reads: 

P>0 => V(p; A) + l -u{p; A) >0, (2.19) 

and then by Lemma 2 we have the possible further consequence : 

detw(p; A) ^ V(p; A) > 0. (2.20) 

Remark: In case the operators A a commute pairwise, in any state p there is a 'classical' 
joint probability distribution over the sets of simultaneous eigenvalues of all the A a . In such 
a case, the term u in Eqs. ( I2.181I2TT91 vanishes identically, and the uncertainty relation (12. 19ft 



is a 'classical' statement 16]. Therefore in the general case a good name for u ab (p; A) is 



that it is the 'commutator correction' term. 

It is instructive to appreciate that while the original definitions of Q and Q, starting from 
the operator sets A and fl, make it essentially trivial to see that they must be nonnegative, 
the form (12. 19[) of the general uncertainty relation gives prominence to the variance matrix 
V(p; A). In addition, as seen earlier, the matrix Q does not directly deal with fluctuations. 
It is after the use of Lemma 1 that we obtain the matrix Cl involving the fluctuations. 

From Eqs. (I2.81I2.17|) . the effect of a unitary symmetry transformation on the real matrices 
V{p ; A) and u{p ; A) is seen to be : 

p' = U pU^ 1 : V(p' ; A) = R V(p ; A) R T , 

u(p'; A) = Ru(p;A)R T , (2.21) 
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so that the form (12. 19j) of the uncertainty relation is manifestly preserved. 

In later work, when there is no danger of confusion, we sometimes omit the arguments p 
and A in V and u. 



III. THE MULTI MODE SCHRODINGER- ROBERTSON UNCERTAINTY PRIN- 
CIPLE 

As a first example of the general framework we consider briefly the Schrodinger- Robertson 
UP for an n-mode system, which has been extensively discussed elsewhere [la . [l7 ]. 

The basic operators, Cartesian coordinates and momenta, consist of n pairs of canonical 
q and p variables obeying the Heisenberg canonical commutation relations. The operator 
properties and relations are : 

?(o+i)/2, a odd, 
p a /2, a even; 

to, = j 



a = 1,2, 



2n : £a 



6a, 6b 

/3 = block diag ( ia 2 , i&2, ■ • • , mt 2 ) = Inxn ® icr 2 . 

These operators act irreducibly on the system Hilbert space H = L 2 (M. n ). 
We take these £ a as the A a of Eq. (12.31) . so here N = 2n : 



Pi 



V ab = - Tr (p {£, - Tr(/)^), 6 - Tr(p&)} 



Then for any state p, the variance matrix V has elements 

1 



(3.1) 



(3.2) 



(3.3) 



while the antisymmetric matrix u> is just the state-independent numerical 'symplectic metric 
matrix' (3 : 



Uab 



) = frfiab- 



(3.4) 
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The uncertainty relation (I2.19P then becomes the n-mode Schrodinger-Robertson UP : 

h 

p > =>■ V + i-f3 > V > 0), (3.5) 

the second step following from Eq. (12.201) as is nonsingular. 

For n = 1, a single mode, the matrices V and (3 are two-dimensional : 



V 



(Aq) 2 A(q,p) 
A(q,p) (Ap) 2 
(Aq) 2 =((q-(q)) 2 ), (Ap) 2 = ((p - (p)) 2 ) , 

&(q,p) = 7;({q- (q),p-(p)}) ; 



Then ( 13. 5ft simplifies to 



(Ag) 2 A(g,p) + fft 
A(q,p)-ih (Ap) 2 



> 0, 



i.e., det \V + -hfij = (Ag) 2 (Ap) 2 - (A(g,p)) 2 - — > 0, 

ft 2 

i.e., det V > — , (3.7) 

the original Schrodinger-Robertson UP. 

Returning to n modes, the Sp(2n, R) covariance of the Schrodinger-Robertson UP (13. 5p 
takes the following form : If S G Sp(2n, R), i.e., any real 2n x 2n matrix obeying S(3S T = (3, 
then the new operators 

ia = S ab ^ b (3.8) 

preserve the commutation relations in Eq. (13. 1 p and hence are unitarily related to the £ Q . 
These unitary transformations constitute the double valued metaplectic unitary representa- 



tion of Sp(2n, R) [18| : 



S E Sp(2n, R) — > U(S) = unitary operator on H, 
U(S')U(S) = ±U(S'S); 

Uffl^LUiS) = S ab ib. (3.9) 
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Then, as an instance of Eqs. ( I2.2ip we have the results : 



p p' = U(S) p U(S)- 1 => V -> V = S V S T , 

V + -h/3>0 V' + -h(3>0. (3.10) 
2 2 



Remark: The n-mode Schrodinger- Robertson UP (13.51) . with its explicit Sp(2n, R) co 



variance (I3.10p . constitutes the answer to an important question raised by Littlejohn [191]: 
under what conditions is a real normalized Gaussian function on a 2n-dimensional phase 
space the Wigner distribution for some quantum state? The answer is stated in terms of the 
variance matrix which of course determines the Gaussian up to phase space displacements 
[And these phase space displacements have no role to play on the 'Wigner quality' of a phase 
space distribution] . This result has been used extensively in both classical and quantum op- 
tics jl^] , and more recently in quantum information theory of continuous variable canonical 



systems 



20|. 



As a last comment we mention that as according to Eq. (13.51) the variance matrix V 



is always positive definite, by Williamson's celebrated t 



found such that V in Eq. (13 .101) becomes diagonal \21 
elements of V will not be the eigenvalues of V. 



reorem an S G Sp(2n, R) can be 



221 ] . In general, though, the diagonal 



IV. HIGHER ORDER MOMENTS FOR SINGLE MODE SYSTEM 



We now revert to the n — 1 case of one canonical pair of hermitian operators q and p, 
but consider expectation values of expressions in these operators of order greater than two. 
The relevant Hilbert space is of course H = L 2 (R). As a useful computational tool we work 
with the Wigner distribution description of quantum states, and the associated Weyl rule 
of association of (hermitian) operators with (real) classical phase space functions. 

Given a quantum mechanical state p, the corresponding Wigner distribution is a function 
on the classical two-dimensional phase space : 



W(q,p) 



1 



2nh 



q ~2 q 



P 



q + ~q') e ivq 'l h 



(4.1) 



Thus it is a partial Fourier transform of the position space matrix elements of p. This 
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function is real and normalised to unity, but need not be pointwise nonnegative : 

pi = p^W(q,p)* = W(q,p); 

/*00 P DC 

Trp = l^ / / dqdpW(q,p) = 1. (4.2) 



— oo J — oo 



The operator p and the function W(q,p) determine each other uniquely. The key property 
is that the quantum expectation values of operator exponentials are equal to the classical 



phase space averages of classical exponentials with respect to W(q,p) [23J : 

/oo poo 
/ dqdpW(q,p)e l(eq - Tp) , -oo < 9, r < oo. (4.3) 
-oo J —oo 

By expanding the exponentials and comparing powers of 9 and r we get : 

nOQ /*00 

Tr(p(q n p n ')) = / dqdpW(q,p)q n p n ', 

J —oo J —oo 

^— ^_ (i9\ n (—ir) n ' 

(qn p n'\ = coefficient of ^ f- in e i( ^" T?3) 

n\ n\ 

77,1 77 I , , 

'—— x coefficient of 9 n (-r) n in (9 q - rp) n+n , 



(n + n')\ 

n, ri = 0, 1, 2, ■ •• . (4.4) 

Thus (q n p n ') is an hermitian operator polynomial in q and p associated to the classical real 
monomial q n p n ' . This is the Weyl rule of association indicated by 

(Ipp) = (q n p n ') w , (4.5) 

so Eq. (14.41) appears as 

TrO%> n V) = / / dqdpW(q,p)q n p n '. (4.6) 

We regard the polynomials (q n p n ') w as the basic 'quantum monomials'. By linearity the 
association (14.51) can be extended to general functions on the classical phase space, leading 
to the scheme : 

f(q,p) = real classical function — y F — (f(q,p))w — hermitian operator on H, 

Tr(pF)= ! fdqdpW(q,p)f(q,p). (4.7) 



Remarks: Two useful comments may be made at this point. For any pair of states p, p' 
we have 

Tr(pp') = J J dqdpW(q,p)W'(q,p) > 0. (4.8) 
12 



Based on this, one can see the following : a given real normalised phase space function 
W(q, p) is a Wigner distribution (corresponding to some physical state p) if and only if the 
overlap integral on the right hand side of Eq. (14. 8p is nonnegative for all Wigner distribu- 
tions W'(q,p). Secondly, we refer to the remarks made following Eq. ( I2.20p concerning the 



commutative case 



A a , Ai 



- 0. This happens for instance when A a = f a (q) for all a. In 
that case, only the integral of W(q,p) over p is relevant, and this is known to be the coor- 
dinate space probability density in the state p [24( . In the multi mode case this generalizes 
to the following statement: the result of integrating W(qi,pi, q%,P2, "• , qn,Pn) over any 
(n-dimensional) linear Lagrangian subspace in phase space is always a genuine probability 
distribution (the marginal) over the 'remaining' n phase space variables. This marginal is 
basically the squared modulus, or probability density in the Born sense, of a wavefunction 



24] 





(a b\ 


H 






dj 



on the corresponding 'configuration space ', generalised to the case of a mixed state 

The covariance group of the canonical commutation relation obeyed by q and p is (apart 
from phase space translations) the group Sp(2, R) : 



real 2x2 matrix | S a 2 S T = a 2 , i.e., det S = 1 > . (4.9) 



The actions on q and p by matrices and by the unitary metaplectic representation of Sp(2, R) 
are connected in this manner : 

S G Sp(2, R) U(S) = unitary operator on % ; 
Z=[\ I = (Ow = f q 

W If 

u{sy 1 iu{s) = si (4.io) 

The effect on W(q,p) = W(£) is then given as j3, Q : 

p' = u(s) pTJ(S)- 1 ^ w'd) = w^s- 1 0- (4.11) 

We now introduce a more suggestive notation for the classical monomials q n p n ' and their 
operator counterparts (q n p n ')w- This is taken from the quantum theory of angular momen- 
tum (QTAM) and uses the fact that finite-dimensional nonunitary irreducible real repre- 
sentations of Sp(2, R) are related to the unitary irreducible representations of SU{2) by 
analytic continuation. (Indeed the two sets of generators are related by the unitary Weyl 
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trick). We use 'quantum numbers' j = 0, ~, 1, • • • , m = j, j — 1, • • • , —j as in QTAM and 
define the hermitian monomial basis for operators on % in this way : 



T jm = (q J+m p 3 - m ) w = coefficient of 



j = o, J,i 

For the first few values of j we have 



(j + m)\(j — m)! 
™ = j, j - 1 



J 5 ) 



0j-H»(_ r )i-"» ^ (Qq- T p) 2j , 
,-■■,-3- (4.12) 



(f 



5 
P 



i ^ lm 



p} 

p2 



\ 



J 



( 



\ 



^(g 2 p + qpq + pq 2 ) 



+ pgp + p 2 g) 
P 3 



/ 



2m 



j(q 3 p + q 2 pq + qpq 2 + pq 3 ) 



\{q 2 p 2 + qpqp + <ff> 2 g + pq 2 p + pqpq + p 2 q 2 ^ 
i :. ,;3 _|_ _|_ p2qp _)_ 



V 



Then we have the consequences : 



P 



(4.13) 



Tr(p7} m ) = y J dqdpW(q,p)q j+m p>- m = qi +m pi 
SESp(2,R): U(S)- l f ]m U(S)= K ml(S)f jm , 



(4.14) 



m——j 



The quantum expectation values of the 7} m are phase space moments of W(q,p), denoted for 
convenience with an overhead bar. The matrices K^'(S) constitute the (2j + l)-dimensional 
real nonunitary irreducible representation of Sp(2, R) obtained from the familiar 'spin j' 
unitary irreducible representation of SU(2) by analytic continuation. For j — |, we have 
K^ 1 ' 2 '^) = S. The representation K^>(S) corresponding to j = 1 will be seen to engage 
our sole attention in Section V. 

The noncommutative (but associative) product law for the hermitian monomial operators 
Tj m has an interesting form, being essentially determined by the SU(2) Clebsch-Gordan 
coefficients. This is not surprising, in view of the connection between 577(2) and Sp(2, R) 
representations (in finite dimensions) mentioned above. In fact for these representations 
and in chosen bases, SU(2) and Sp(2, R) share the same Clebsch-Gordan coefficients 
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The product formula has a particularly simple structure if we (momentarily) use suitable 
numerical multiples of Tj m : 

(4.15) 



Tj m = f jm / a/ (j + mjUj - m)\. 



Then we find [14 ] 



Tjm Tj'm' 



j"=\j-f\ 
•jj'f 

m m 



(2j» + l)(j + f - f)\(f + j" - j)\(j" + j - 



y^fj j j _ 

X U , m+m , Tj" : m +m' ■ 



(4.16) 



The C^^, „ are the SU(2) Clebsch-Gordan coefficients familiar from QTAM 25J. We will 
use this product rule in the sequel. 

Now we apply the general framework of Section II to the present situation. We will use 
a notation similar to that in the main theorems of the classical theory of moments. We take 
A and A to formally be infinite component column vectors with hermitian entries : 
/ : \ 



A 



A 



T 



\ ' J 
1 

A 



(f 



ii) Ti=i, Tii, Tio, Ti _i, 

2 2 2 2 



T- ■ 



T ■ 

' ± 3,-3i 



(4.17) 



Thus the subscript a of Eq. (I2.4p is now the pair jm taking values in the sequence given 
above. To simplify notation, as A is kept fixed, we will not indicate it as an argument 
in various quantities. The general entries in the infinite-dimensional matrices Cl, Q, Q in 
Eqs. (El EHD are then : 



%tnj')n' Tj m Tji m i , 



'(p) — Ti(pTj m Tj, n/ ' 



\Tjm Tj'm' ) 



Qjm,j'm'[P) — ((Tjm — (Tjm)) {Tj'm' ~ (Tj'm'))) ■ (4-18) 

(In Cl and Q, for j = m = 0, we have T 00 = 1). By using the product rule (14.161) the 



(generally nonhermitian) operator Tj m T,y m t can be written as a complex linear combination 
of Tj'^m+m' with j" = j + j', j + j' — 1, • • • , \j — f\. The variance matrix V(p) in Eq. ( 12.181) 
has the elements 



*jm,j'm' 



(P) = -^({Tjm, Tj' m '}) - (Tj m ) (Tj' m >) . 



(4.19) 
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From the known symmetry relation 25] 



rjj'ii" _ (_i)j+3'-i" nii'i" (a oq) 

in' mm+m 1 V / mm' m+m' V " w / 

we see that in the anticommutator term in Eq. (14.190 only T m+m , for j" = j + f, j + j' — 
2, j — 4, • • • will appear with real coefficients. On the other hand, for the antisymmetric 
part u a b of Eq. (12.180 we have 

Ujm,j>m>(f>) = — i ( Tj m , Ty m i ) , (4-21) 



so now by Eq. (14.200 the commutator here is a linear combination of terms T^ +m , for j" = 
j + j' — 1, j + j' — 3, • • • with pure imaginary coefficients. There is, therefore, a clean 
separation of the product Tj m Tj/ m i into a hermitian part in V and an antihermitian part in 
oj. With these facts in mind, the uncertainty relation (12.190 is in hand: 

j+j'—j" even 

(P) = " ' ( T f',rn+m') ', 

j+j'-j" odd 



bJjm,j'm' 



(%m,i'm'(p)) = (Vjm,j'm'(p)) + - (^mj'm'(p)) > 0. (4.22) 

Each matrix element of V(p) (apart from the subtracted term) and of oo(p) appears as 
some real linear combination of expectation values of hermitian monomial operators, i.e., of 
moments of W(q,p) ; however, in this way of writing, the essentially trivial nature of the 
statement fl(p) > is not manifest. 

The covariance group in this problem is of course Sp(2, R). From Eq. (14.140 we see that 
under conjugation by the metaplectic group unitary operator U(S), the column vector A of 
Eq. ( 14.170 transforms as a direct sum of the sequence of finite-dimensional real irreducible 
nonunitary representation matrices K^^ 2 \S) = S, K^\S), K^ 3 ^ 2 \S) ■■■ ; so Eq. (12.80 in 
the present context is Jjj] : 

S e Sp{2, R) : U(S)- 1 AU(S) = K{S) i, 

K{S) = K il/2 \S) © K {1) (S) © if (3/2) (S) © • • • (4.23) 

From Eq. (12.211) . when p — > p' = U(S) pUlS) -1 both V(p) and u(p) experience congruence 
transformations by K(S), and the formal uncertainty relation (14.221) is preserved. 
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Up to this point the use of infinite component A and infinite-dimensional Q, fi, V and 
to has been formal. We may now interpret the uncertainty relation f)4.22p in practical terms 
to mean that for each finite iV = 1, 2, • • • , the principal submatrix of Cl(pi) formed by its 
first N rows and columns should be nonnegative. However, in order to maintain Sp(2, R) 
covariance, a slight modification of this procedure is desirable. If for each J = ~, 1, |, • • • 
we include all values of j m for j < J, the number of rows (and columns) involved is 
Nj = J[2J + 3), the sequence of integers 2, 5, 9, 14, ■ ■ • . Let us then define hierarchies of 
iVj-dimensional matrices as : 

J= r 1, 2'"' : 

^ (J) (P) = (^mj'm'(p)), 
V (J \f>) = (V jm ,fm>(p)), 

i 

U (J \p) = {Ujm,j>m'{p)), j, 3 

^ J \p) = V^ J \p)+ l -^ J \p). 

However, in each of these matrices there is no J dependence in their matrix elements. Each 
also naturally breaks up into blocks of dimension (2j + l)(2j' + l) for each pair (j, j') present, 
and these can be denoted by Cl^(p), V^(p), u^ j '\p). Symbolically, 

/ \ 

& J \p) = 



2' 



(4.24) 



\ 



(4.25) 



/ 



Q^(p) 



(4.26) 



and similary for V^ J \p) and cu^ J \p). As examples we have: 

n^(p) = (n(^)(p)); 

(p) n^ip) 
n(^)(p) nW(p) 

and correspondingly for V^ J \ to^ J \ Moreover, in going from J to J + |, we have an aug- 
mentation of each matrix with 2( J + 1) new rows and columns, 

/ ; n(i.*4)(p) \ 

fi (J) (p) ; ; 

^(J+i/2) (/3)= : a(J,J+h)(fi . (4.27) 



] (p) 
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The formal uncertainty relation (I4.22p now translates into a hierarchy of finite-dimensional 
matrix conditions 

& J \p) = V^{p) + l -^ J \p) >0, J=i 1, ? (4.28) 

(Of course, for a given state p, moments may exist and be finite only up to some value J max 
of J, so the hierarchy (14. 28ft also terminates at this point). The lowest condition in this 
hierarchy, J = |, takes us back to Eqs. (13.61 13. 71) : 



2 



(Aq) 2 A(q,p) 
A(q,p) (Ap) 2 ^ 

u^\p) = Lo^'^(p) = -i ( ^ ^' ^ | =iha 2 

\\P,q] o y 
ft(i)(fl>o * f (A?)2 A( ^ ) + ^ 

A(q,p)-il (Ap) 2 



> <£> 



(Ag) 2 (Ap) 2 -(A(g,p)) 2 >^, (4.29) 

the original Schrodinger-Robertson UP. 

It is natural to ask for the new conditions that appear at each step in the hierarchy (14.281) . 
in passing from J to J + |. In the generic case, when we have a strict inequality we can find 
the answer using Lemma 1 of Section II. Comparing Cl^ J+ ^\p) and ^ J \p), in the notation 
of Eq. (l2~T0j) and using Eq. <K27} we have : 



fi( J+ %) 



A C f 
C B 

A = & J \p), B = n( J+1 i' J+ ^)(p), 

C = ( n( J+U) ( p ) ... c^+M (p) ) . (4.30) 
The 'dimensions' are Nj x iVj, 2( J + 1) x 2( J + 1), 2( J + l)xiVj respectively. Then 

fi( J+ ^)(p) > fi (J) (^) > 0, B-CA- l Cfl>0, (4.31) 

18 



where A, B, C are taken from Eq. (14.301) . One can see that some complication arises from 
the need to compute A~ l in the new condition. 

In the next Section, we analyse the case J=|— 7-J+| = lin some detail, as the first 
nontrivial step going beyond the Schrodinger- Robertson UP (13.71 14.291 . Before we turn to 
this task, however, a note on the non-generic case of singular A seems to be in order. 

Remark : Lemma 1 expresses the positive definiteness of a hermitian matrix Q in the block 
form (12.101) in terms of conditions on the lower dimensional blocks. The block form itself is a 
description of Q with respect to a given breakup of the underlying vector space on which Q 
acts, into two mutually orthogonal subspaces. Both A and B are hermitian. For the case of 
positive semidefmite Q, there are two possibilities at the level of A, B, C. If A" 1 exists, then 
Q > translates into A > 0, B - C A' 1 C f > 0. In case A is singular, while of course Q > 
implies A > 0, the question is what other condition on B, C is implied. To answer this, we 
further separate the subspace on which A acts into two mutually orthogonal subspaces — one 
corresponding to the null subspace of A, and the other on which A acts invertibly, say as 
A\. Then in such a description, the block form of Q is initially refined to the form 

^ C\ ^ 

Q- At C\ , (4.32) 

yC 2 Ci B J 

with the original A and C being respectively and (C2, CM. But now one sees 

V' A J 

easily that Q > implies C 2 = 0, so as A 1 1 exists, we have in this situation 

Q > Ax>0, B-dA^ 1 C\ > 0. (4.33) 
This is the description of the nongeneric situation mentioned above. 



V. SO(2, 1) ANALYSIS OF FOURTH ORDER MOMENTS 

The first nontrivial step in the hierarchy of uncertainty relations (I4.28p . after the 
Schrodinger- Robertson UP (13. 7\ 14.291) . occurs in going from J=|toJ+^ = l. We study 
this in some detail, especially as it brings into evidence the equivalence of the irreducible 
representation K^(S) of Sp(2, R) and the defining representation of the three-dimensional 
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proper homogeneous Lorentz group SO(2, 1) |26j. Indeed K^ 2 >(S), K( 3 \S), ■ ■ ■ are all true 
representations of 5*0(2, 1) |27|. 

It is useful to introduce specific symbols for the operators Ti m , T lm in the present context. 
We write 



{Ti m ) = (im) 



{Tim) — {X r , 




m — -, — ; 
2' 2 ' 



5 {?> P} 



m — 1, 0, — 1 ; 



(5.1) 



/ 



so that one immediately recognises that £ is a two-component spinor, and X a three- 
component vector, with respect to 50(2, 1) (see below). Their products can be computed 
using Eq. (14.161) or more directly by simple algebra : 



irnim' — X rn+m i + % | ( — l) m 2 8, 



m+m' , j 



£m X n 



2 



, +i 5(_l)m- 3 £ 



h 

, m+m' 1 " 2 
/I 2 



m+m 1 i 



(b) 



(5.2) 



X m X m / = T 2 ,m+ m > + -H-l) m (l +m 2 )5m+ m ',o + i h (m - m') X m+m , . (c) 



In (5.2a) the leading J = 1 term is symmetric in m, m' ; while the pure imaginary J = 
second term is antisymmetric. In (5.26) it is understood that £ ± 3 = 0. In (5.2c) the first two 
J = 2 and J = terms are symmetric in m, m' ; while the third J = 1 term is antisymmetric. 
These features agree with the pattern in Eq. f)4.22p . 
For J = | in Eq. (I4.30P we have 



(5.3) 



with 'dimensions' 2x2, 3x3, 3x2 respectively. (Throughout this Section, A, B, C will 
have these meanings). Their behaviours under Sp(2, R) are 

p^p' = U(S)pU(S)- 1 A^SAS T , B KW(S)BKW(S) t , 

C K (1 \S)C S T . (5.4) 



Assuming A 1 exists, we have 



A- 1 ^ (S- 1 ) 7, A- 1 S~\ 



(5.5) 
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and consequently, 

B - C A' 1 O f -> K {l \S) (B -OA' 1 O f ) if (1) (£) T , 



(5.6) 



which as expected is a congruence. 

The matrix K^(S) is easily found. At the level of classical variables: 



$ = I a b I G £p(2, i?) : I 9 I -)• S 1 H 
c d } \ p I \ p 



{X m (q, p)) 



-^«(£) (X m (g,p)), 



( a 2 lab b 2 ^ 



K (1) (£) 



V 



ac ad + 6c 6d 
2cd d 2 



(5.7) 



/ 



The link to 50(2, 1) can be seen in two (essentially equivalent) ways, either through A or 
through (X m (q,p)). We now outline both. 

We introduce indices fi, u, ■ ■ ■ going over values 0, 3, 1 (in that sequence) and a three- 
dimensional Lorentz metric g^ v = diag(+l, —1, —1). This metric and its inverse g^ v are 
used for lowering and raising Greek indices. The defining representation of the proper 
homogeneous Lorentz group £0(2, 1) is then: 

£0(2, 1) = {A = (A"„) = 3x3 real matrix | A M = 9flT A"„ A T A = g uX , 

detA = +1, A% > 1}. (5.8) 

This is a three-parameter noncompact Lie group. Now expand A = (p) in terms of 

Pauli matrices as follows : 



X° + X 3 x 1 



h 



X 



-cr 2 - 



(5.9) 



From Eqs. (13.61 14. 29j) we have (indicating p dependences) : 

x°(p) = l -{{Aqf + (Ap) 2 ), x 3 {p) = \{{Aq) 2 - (Ap) 2 ), x\p) = A(q,p) (5.10) 
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Then the transformation rule for A in Eq. (15.40 , combined with S 02 S T = <Ji , leads to a rule 
for x^{p) : 



P^U(S)pU(S)- 1 =4> A^SAS T => x»(p) A» u (S)x u (p), 



( 1 f a 2 + & 2 + c 2 + d 2) 1 (a 2 _ 6 2 + c 2 _ d 2N a& + crf \ 



A(S) 



V 



±(a 2 + fe 2 -c 2 -d 2 ) i( 
ac + fed 



ac — fed 



+ d 2 ) afe — cd 
ad + fee 



G 50(2,1). (5.11) 



Thus x^{p) transforms as a Lorentz three- vector, and the associated invariant is seen to be 

x^p)x,{p) = 9llv x»(p)x»(p) = (Aq) 2 (Ap) 2 - (A(q,p)) 2 > £ (5.12) 

so the Schrodinger- Robertson UP implies the geometrical statement that x^(p) is positive 
time-like. 

The matrices K^(S) by which X m transform under Sp(2, R) are related by a fixed 
similarity transform to the A(S) above. If in terms of classical variables we pass from the 
components X m (q,p) in Eq. (15. 7ft to a new set of components X^(q,p) by 









(A 




W-p 2 ) 


= M 


qp 




V qp ) 




\ P2 J 



X»(q,p) = M» m X m (q,p), X m (q,p) = M^X»(q,p), 



/1 n i\ 



M = (M p 



2 2 

2 u 2 
1 



M- 1 = (M-J) 



(1 1 o^ 

1 

1 -1 



then in place of Eq. ( 15.70 we have 



(5.13) 



Q 
P 



q 
p 



X"{q,p) M^ m K { J l l(S)M-} u X%q,p) 

= A^(S)X»(q,p), 
K {1) {S) = M~ l A(S) M. 



(5.14) 



At the operator level we have 

1 



i(g 2 -p 2 ), X^ X -{q,p}, 



(5.15) 
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and, as consequence of Eq. (I4.14p ). the twin equivalent transformation laws: 



S e Sp(2, R) : U(S)- 1 X m U(S) = K { ±,(S)X, 
U(S)- 1 X^U(S) = A^S) X v . 



(5.16) 



The upshot is that the matrices K^(S) are just the 'ordinary' homogeneous Lorentz trans- 
formation matrices A(S) in a 'tilted' basis. The metric preserved by them is easily found 
though unfamiliar : 



K^\S)g K K^\S) T = g K) 



( n n o \ 



g K = M- 1 g(M- 1 ) 



1\T 



2 
0-10 
2 



(5.17) 



This enables us to use the nomenclature and geometrical features of three-dimensional 
Minkowski space even while working with operators X m and transformation matrices 
K&(S). 

Now we proceed to analyse the three matrices A, B, C and the combination B — C A^ 1 C'. 
(We have already parametrised A in Eqs. (15 .9[ 15.101) ). Using Eqs. (I5.2p . their matrix elements 



are 



.4, 



(Cm Crn') (Cm) (Cm') 

(x m+m ,) - (u d m ,) + 1 ? (-i) m -^ s r , 



h 



m, — m j 



B r , 



C,, 



(X m X m r) — (X m ) (X m r) 

h 2 ... . 

(T2,m+m>) + —T ( — l) m 0~m,-m' ~ (X m ) (X m r) + % fl [m — m ) (X m+m f)] 



A 



\-X m Cm') (X m } (Cm') 



h 



( T f ,m+m') _ ( X m) (Cm') ~ % ^ 2 (Cm+m')- 



(5.18) 



In each of these expressions, the possible values for m, m' are evident from the context. 
We now note an important fact in respect of the final forms of all three expressions: apart 
from explicit appearances of i in the last terms, all other quantities are real. This allows us 
to easily separate each of A, B, C into real and imaginary parts, which in the cases of A 
and B are respectively symmetric and antisymmetric in m and m'. [This is already seen in 
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Eq. (15. 9 p for A}. We write these as follows : 

A = A 1 +iA 2} 

(-^Omm' V^m+m') (Cm) (Cm'/ (x^ 4 0"^) mm ' , 

B = B 1 +iB 2 , 

(^l)mm' = (^2,m+m') + IT ( — 1)™ $m,-m' ~~ (X m ) (X m i) , 

(5.19) 

(■B 2 )mm' = h{m-m') (X m+m >); (b) 
C = Ca+iC 2 ; 

(Cl)mm' = (7| )m+m ') - (^m) (Cm'}, 
(C 2 )mm' — ~ | ( — 1)™ 5 (Cm+m') j 

Ct = Cf-iCj. (c) 

To deal similarly with B — C A^ 1 C\ we need an expression for A -1 . We will assume the 
generic situation in which A is nonsingular, 

h 2 

detA = K- 1 =x^x fl - — >0, (5.20) 



so that 



h 

A' 1 = k(x° - x 3 (T 3 - iVi + - tr 2 ) 
= K^arap + - a 2 ), 



(x°, -x 3 , 



-x 1 ). (5.21) 



The transformation law for A 1 under S G Sp(2, R) given in Eq. (15. 5 j) is different from 
(though equivalent to) the law for A. Thus, while the x M do follow a definite (i.e., well 



defined tensorial [26| ) transformation law, there are some differences (in signs) compared 
to the law followed by x^ 1 . Clearly the two terms in Eq. (15.211) are, as they stand, the 
real symmetric and the pure imaginary antisymmetric parts of A~ l . We can now handle 
B — C A^ 1 in the same manner as above : 

B - C A' 1 C ] = B 1 + % B 2 - «(d + iC 2 )(x-a + ^ a 2 ) (Cf -%C T 2 ) 

= y(eff) + i u {eff) 

2 

= B 1 -K(C 1 x-aC? + C 2 x-*C 2 T +i^C 2 a 2 Ci;-i^C 1 a 2 C 2 T ), 
l - u {cS) =B 2 - k(C 2 x-aCl-dx. aCj - % -C x a 2 C T X - i- C 2 o 2 C T 2 ). (5.22) 
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This decomposition is in the spirit and notation of Eq. (I2.18P of the general framework. 
However, I/( eff ) and u/ cff ) do not correspond any longer to expectation values of sim- 
ple anticommutators and commutators among relevant operators, as was the case in 
Eqs. fl2M8[l4.i9[IPT|) . 
Both and 

^(cff) are rea j three-dimensional matrices with elements , oo^^i, where 
m, m' = 1, 0, —1; and they are respectively symmetric and antisymmetric. It does not 
seem possible to simplify the expressions (15. 22ft to any significant extent, as they are already 
expressed in terms of the independent real expectation values (£ m ), (X m ), (T| m ), {T 2 ^ m ) 
which are the moments of the Wigner distribution W(q, p) of orders up to and including 
the fourth. Under action by S G Sp(2, R) we have from Eq. (15. 6 P : 

p^U(S)pU(S)- 1 V {cS) ^ K {1 \S)V {cS) K^iSf, 

^K^(S)u^K^(Sf. (5.23) 

The added uncertainty relation up to the fourth order going beyond the Schrodinger- 
Robertson UP (13.71 14.291) . reads [in the generic case det A > 0] : 

y(eS) + % ^(eff) > Q) (534) 

which is an SO(2, 1) covariant statement by virtue of Eq. (15.231) . 

For further analysis it is rather awkward to work with SO(2, 1) matrices and Lorentz 
metric in the form K^(S), gx, therefore we pass to the 'standard' forms via the matrices 
M, M" 1 in Eq. (1533]) : 

VWi» = M» m M v m ,v£S, 

^'"^M^MVwS. (5-25) 
which are congruences. Then the Sp(2, R) or SO(2, 1) actions (I5.23P appear as: 

y(eS)^ _^ AV(^) A" AS) U (eff)MV ', 

jcs)^ _^ AVOS) A\,{S) w (eff ^ v , (5.26) 
and the condition (I5.24p becomes : 

(y( eff )^) + i(a;( efr )^) > o. (5.27) 
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While ^/( eff )^ transforms as a symmetric second rank 5*0(2,1) tensor, anti- 
symmetric second rank tensor, which by the use of the Levi Civita invariant tensor is the 
same as a three vector. Thus we can write, with e 031 = e 3i = +1, 



with transformation law 



(CO 



(eff)ju/\ 



a% —a 3 * 
— a\ ao 
03 — ao 



(5.28) 



a" -> A^(5) a 1 '. 



(5.29) 



Of course, ^ v itself is made up of two irreducible parts : the symmetric second rank 
'trace-free' part belonging to the 50 (2, 1) representation K^(S), and the 5*0(2, 1) invariant 
trace which is a scalar. 

We now appeal to a remarkable result j^ . which is similar in spirit to the Williamson 
theorem for Sp(2n, R) quoted in Section 3. It states that if ^/( eff )^ transforming as in 
Eq. (I5.26P is positive definite, it is possible to bring it to a diagonal form by a suitable choice 
of A G 50(2, 1) ; however in general the resulting diagonal values are not the eigenvalues of 
the initial matrix. This diagonal form may be called the 'SCS normal form' of V^ eS \ which 
in the generic case is unique. Passing to this normal form of V^ eS \ and transforming u/ eff ) 
as well by the same (generically unique) A G 50(2, 1), these matrices appear as 



y(eS) 



— >■ 





v 33 
v 11 



(eff) 



— > 



\ 



-6 1 b 3 
b 1 6° 
-b 3 -b° 



(5.30) 



with all the quantities v 



00 ^33 ^11 



6°, b 3 , b 1 being real 50(2, 1) (and Sp(2, R)) invariants. 



The uncertainty relation f !5.27p expressed in terms of these invariants, and in its maximally 
simplified form thanks to the SCS theorem, is 





v 33 
v 11 



\ 



/ 



+ 






-b 1 


b 3 


b 1 





b° 


-b 3 


-6° 






> 0. 



(5.31) 
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As an (admittedly elementary) example of the discussion of this Section, we consider the 
Fock states |n), n > 0. The (q, p) — (a, a^) relations are 

a = (q + ip) /V2h, a) = (q- ip) /V2h, (5.32) 

so both q and p have dimensions h 1 ^ 2 . In the Fock states |n), by parity arguments we have 

(n\i m \n) = (n\f 3/2 ,m\n) = 0. (5.33) 

For X m , T 2 ,m explicit calculations give : 



(n\X m \n) = h(n + -)(1, 0, 1), m = 1, 0, -1; 
(n\f 2tm \n) = h?{n 2 + n + ~)(3, 0, 1, 0, 3), m = 2, 1, 0, -1, -2. 



Then the matrices A, B, C of Eq. ( 15.31) follow easily: 



(A 



/ 



2> 2 



- % + i), x 3 = x 1 = 0; 



f (n 2 + n + 1) 



0. 



1 -1 \ 
1 
-10 1 / 



+ i^ 2 (n + i; 



/ 



V 






1 





-1 





1 





-1 






(5.34) 



(b) 



(5.35) 



Therefore the combination B — CA = B, and from Eq. (15.22ft . 



mm 



h 2 



[n 2 + n + 1) 



V 



2 V mm 



h 2 (n + -) 



1 
-1 1 
0-10 



1 -1 \ 

1 

1 1 

\ 



(5.36) 



Transforming to the standard 50(2, 1) tensor components by the congruence transforma- 
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tion ( I5.25P we find : 



^0 ^ 



(n 2 + n + 1) 



h 2 (n+-) 



1 

1 

0^ 

1 

1 0, 



(5.37) 



As expected, both these matrices are invariant under the 5*0(2) subgroup of 5*0(2, 1), as 
the Fock states are eigenstates of the phase space rotation generator a^a. 

We see that (y^^) is already in the SCS normal form, and as the eigenvalues of 



(yW^+i (o/eff)/^) ar e 0, \(n 2 + n + 1) ±h 2 (n + §), i.e., 0, \{n + l){n + 2), \n{n-l), 
the uncertainty relation (I5.27P is clearly respected; indeed it is saturated! 



VI. LORENTZ GEOMETRY AND THE SCHRODINGER- ROBERTSON UP 

The original Schrodinger-Robertson UP has a very interesting character when viewed in 
the Wigner distribution language, bringing out the role of the group 50(2, 1) in a rather 
striking manner. This seems worth exploring in some detail. 

For a given state p with Wigner distribution W(q, p), the means are 




dqdpqW(q, p) , p 




dqdppW(q, p) . 



(6.1) 



Referring to Eq. (|5.13p . at each point (q, p) in the phase plane we define the 50(2, 1) three- 
vector (a displaced form of (X^(q,p)) in Eq. ( I5.13P ) 



X»(q, p) 



\ 



\[{q-q? + (v-v?} 

\ (?-?)(p-P) j 



(6.2) 



which (except at q = q, p = p) is pointwise positive light-like. The elements of the variance 
matrix V in Eqs. (3.6, 4.29) are obtained by 'averaging' this three-vector over the phase 
plane with the quasiprobability W(q,p) as 'weight' function, resulting in the three- vector 
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x»(p) of Eq. f l5U0|) 




Given that W(q,p) can in principle be negative over certain regions of the phase space, 
this 'averaging' could have led to a result which need not be either time-like or light-like 
positive. However the Schrodinger-Robertson UP assures us that in fact the result has to 
be a time-like positive three-vector, thus implying a subtle limit on the extent to which 
W(q,p) could become negative. In fact it specifies that the three- vector obtained as a 
result of the 'averaging' must be within or on the positive time-like (solid) hyperboloid 
^2 x M (p)x At (p) > h 2 /4 corresponding to 'squared mass' h 2 /4 presented in Eq. (I5.12p . On the 
other hand, while pointwise nonnegativity of W(q,p) will certainly ensure that the averaging 
in Eq. ( 16. 3 p takes (x M (p)) inside the time-like positive cone, it will not itself ensure that it is 
taken all the way inside the said hyperboloid. To ensure the latter, W(q,p) should have a 
threshold effective spread. Thus, pointwise nonnegativity is neither a necessary nor sufficient 
requirement to ensure 'Wigner quality' on W(q,p) as is known from other considerations. 

The argument above has been presented in such a way that the interpretation in terms of 
Lorentz geometry in 2 + 1 dimensions is obvious. However, comparing Eqs. ( 15.71) and (15.131) . 
we see that it could be expressed equally well as follows. At each point (q,p) in the phase 
plane we define a 2 x 2 real symmetric matrix 



Pointwise (except at q = q, p = p ) this is proportional to a one-dimensional projection 
matrix, and in particular it has vanishing determinant. After 'averaging' with W(q,p) as 
weight function, however, we obtain the 2x2 variance matrix V in Eq. ( I4.29P : 



and now the Schrodinger-Robertson UP shows that V is non-singular and has determinant 
bounded below by the 'squared mass' h 2 /4. 

In this form, just like the Schrodinger-Robertson UP, this geometrical picture based on 
the Wigner distribution language generalises in both directions — second order moments for 




(6.4) 




(6.5) 
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a multi mode system, and higher order moments for a single mode system. As an example 
of the former, consider a two-mode system for simplicity. The classical phase space variables 
are £ a and the hermitian quantum operators obeying Eq. (13.1 j) are £ a , for a = 1, • • • ,4. Given 
a two-mode state p, we pass to its Wigner distribution W(£) (something we did not do in 
Section III) and compute the means 

(O = Tr(pi a ) = J d 4 ^ a W(0 = l„ a = 1, • • ■ ,4. (6.6) 

Then, generalising Eq. (16. 4 j) above, at each point £ in the 4-dimensional phase space we 
define a real symmetric 4x4 matrix 

f(0 = (K&(0) = ((& - U(6 - £&)) = *(0*(e) T , 

At each point £ (except at £ = £) we have here a real symmetric positive semidefinite matrix 
V(£) which is essentially a one-dimensional projection matrix: the eigenvalues of V(£) are 
x(£) T x(£), 0, 0, 0. The variance matrix V for the state p is then obtained by 'averaging' V(£) 
using the real normalised quasiprobability : 

v = | A^(£M0- (6-8) 

Since in general W(£) can assume negative values at some places in phase space, it may 
appear at first sight that some of the properties of V(£) described above may be lost by 
the 'averaging' process leading to V. However the UP ( 13.51) guarantees that this will not 
happen; indeed by Lemma 2, Section II, in Eq. (I2.13p . V is seen to be positive definite. 
Quantitatively we have the following situation : Williamson's theorem assures us that under 
the congruence transformation by a suitable So £ Sp(4,7V), V is taken to a diagonal form: 

Vq = SqVSq = diag(ni, k 1: k 2 , k 2 ), k 1)2 > 0. (6.9) 



The congruence transformation becomes a similarity transformation onV(3 1 13J, since: 

S G Sp{4, K) : V' = SVS T ^ V'^ 1 = SV^S' 1 . (6.10) 
Applying this to the transition V — > Vq, we see that as 

, K\<J 2 . 

VoP- 1 = ~i\ | , (6.11) 

K 2 a 2 
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the eigenvalues of iVfl 1 are ±ki, ±^2; and the UP (13.51) ensures that K\ } 2 > h/2. The k's 
themselves are determined, upto an interchange, by the Sp(A, TV) invariant traces 

Tr(y/3- 1 ) 2 = -2(«? + ^), 
TriVr 1 ) 4 = 2(4 + 4). (6-12) 

The manner in which the geometrical picture, and the constraint on the extent to which 
W(£) can become negative, both generalise in going to the multi mode situation is now 
clear. 

A qualitatively similar situation (even if algebraically more involved) obtains when we 
generalise in the other direction — to higher order moments for a single mode system, and 
their uncertainty relations handled in the Wigner distribution language. Limiting ourselves 
to the moments upto fourth order, we are concerned in the notation of Eq. (I4.24p with the 
uncertainty relation 

Q (1) (p)>0 (6.13) 

contained in the hierarchy (14 .28 p . and its rendering in the Wigner distribution language. 
Combining the notations of Sections II and V, we have a set of five real phase space functions 
A a (q,p), a = 1, 2, • • • ,5, and their hermitian operator counterparts in the Weyl sense : 

(A a (q,p)) = (q, p, q 2 , qp, p 2 ) T ; 



(A a ) = ((A a (q,p)) w ) = (q, p, q , -{q,p}, p ) 



2\T 



(6.14) 



a listing of the components | m , X m . In a given state p with Wigner distribution W(q,p) we 
have the means 

(A a ) = Tr(pi a ) = J J dpdqW(q,p)A a (q,p) = A a . (6.15) 

To calculate the elements of Q^(p) we need to deal with the products A a Ab. For these, 
using Eq. ( 15.2ft we find : 

A a A b = (A a (q,p)A b (q,p)) w + {C ab (q,p)) w , 



(Cab(q,p)) 



/ 
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ihp 
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(6.16) 
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(We note that the real symmetric part of the matrix C(q,p) is — \qk in the lower 3x3 
block, where gx is the tilted form of the (2 + 1) Lorentz metric in Eq. (I5.17P ). With these 
ingredients and referring to the general structure (I2.16P we have the expression for Cl^(p) 
in the Wigner distribution language : 

= (n£?(p)) = (Tr(p(i a - (A a ))(A b - (i 6 )))) 

= (Tr(/5 ( (A a (q,p)A b (q,p)) w - A a A b + {C ab (q,p)) w ) )) 

= J J dpdqW(q,p) (x(q,p)x(q,p) T + {C ab (q,p))) , 

x {QiP) T — (<? — Qi v — Pi i 2 — q 2 : qp — p 2 ~ p 2 )- (6-17) 

At each point (q,p) in the phase plane, we have essentially a one-dimensional projector 
x(q,p)x(q,p) T , together with a five- dimensional hermitian matrix {C ab (q,p)) with elements 
involving fi and h 2 terms. The uncertainty relation ( 16. 131) demands that the phase plane 
'average' of this expression (hermitian matrix) with W(q,p) as weight function be nonneg- 
ative. After this 'averaging', the leading term is no longer a one-dimensional projector; 
moreover, the pure imaginary antisymmetric part coming from this part of C(q,p) being 
singular, Lemma 2 of Section II does not apply to the real symmetric part of Cl^(p). In any 
event, (I6.13P again constrains the extent to which W(q,p) can become negative. 

VII. CONCLUDING REMARKS 

In this paper we have set up a systematic procedure to obtain covariant uncertainty 
relations for general quantum systems. It applies equally well to continuous variable systems 
and to systems described by finite-dimensional Hilbert spaces, and even to systems based on 
the tensor product of the two, and consists of two ingredients : the choice of a collection of 
observables, and the action of unitary symmetry operations on them. We have shown that 
the uncertainty relations are automatically covariant — preserved in content — under every 
symmetry operation. 

We have applied this to two important special cases: the fluctuations and covariances 
in coordinates and momenta of an n-mode canonical system; and to the set of all hermi- 
tian operator 'monomials' in canonical variables q, p of a single mode system. These are 
both generalisations of the Schrodinger-Robertson UP in two distinct directions. The latter 
generalisation has been treated for definiteness using the Wigner distribution method. 
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We hope to have set up a robust yet flexible formalism which can be applied to all 
quantum systems, in particular to composite, for instance bipartite, systems. In such a 
case, by judicious choices of the operator sets {A a } of Section II, one can devise tests for 
entanglement, exhibiting covariance under corresponding local symmetry operations. If for 



a bipartite system the operator p PT [28|, arising from partial transpose of a ph ysic al state 



p, violates any uncertainty relation, the presence of entanglement in p follows [29l. [30j. A 
systematic analysis along these lines of higher order moments in the bipartite multi-mode 
case will be presented elsewhere, keeping in mind that our general methods are applicable 
for both discrete and continuous variable systems, and even to composite systems consisting 
of either or both types as subsystems. 
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